Exploration of physiological and biochemical processes of canola with exogenously applied fertilizers and plant growth regulators under drought stress

Environmental stresses may alter the nutritional profile and economic value of crops. Chemical fertilizers and phytohormones are major sources which can enhance the canola production under stressful conditions. Physio-biochemical responses of canola altered remarkably with the use of nitrogen/phosphorus/potassium (N/P/K) fertilizers and plant growth regulators (PGRs) under drought stress. The major aim of current study was to evaluate nutritional quality and physio-biochemical modulation in canola (Brassica napus L.) from early growth to seed stage with NPK and PGRs in different water regimes. To monitor biochemical and physiological processes in canola, two season field experiment was conducted as spilt plot under randomized complete block design (RCBD) with four treatments (Control, Chemical fertilizers [N (90 kg/ha), P and K (45 kg ha-1)], PGRs; indole acetic acid (IAA) 15g ha-1, gibberellic acid (GA3) 15g ha-1 and the combination of NPK and PGRs] under different irrigations regimes (60, 100, 120, 150 mm evaporations). Water stress enhanced peroxidase (POD), catalase (CAT), superoxide dismutase (SOD), polyphenol oxidase (PPO), soluble sugar, malondialdehyde (MDA), proline contents as well as leaf temperature while substantially reduced leaf water contents (21%), stomatal conductance (50%), chlorophyll contents (10–67%), membrane stability index (24%) and grain yield (30%) of canola. However, the combined application of NPK and PGR further increased the enzymatic antioxidant pool, soluble sugars, along with recovery of leaf water contents, chlorophyll contents, stomatal conductance and membrane stability index but decreased the proline contents and leaf temperature at different rate of evaporation. There is positive interaction of applied elicitors to the water stress in canola except leaf area. The outcomes depicted that the combination of NPK with PGRs improved the various morpho-physiological as well as biochemical parameters and reduced the pressure of chemical fertilizers cost about 60%. It had also reduced the deleterious effect of water limitation on the physiology and grain yield and oil contents of canola in field experiments.


Introduction
Among oilseed crops, canola (Brassica napus L.) holds a special rank worldwide as the second highly significant oil yielding crop. with over~40-50% oil and 40% protein (in rapeseed meal) contents after soybean [1]. Owing to continuous efforts of scientists, the production of canola had reached a milestone by enhancing its yield and area upto 68.9 million metric tons and 33.7 million hectares, respectively [2]. The leaves of this plant are imperative source of animal feed due to encompassing a balanced ratio of protein and fibrous [3]. Oil extracted from oilseed rape contains a very high concentration of unsaturated fatty acids (such as oleic acid and linoleic acid (C18:3,~10% v/v)) making it a fatty acid-rich diet source [4]. Pakistan facing the low yield of canola somewhat due to the poor management of macro-and micronutrient. Globally, it is emerging oil seed crop among the other oil seeds due to low erucic acid and glucosinolates in oil and seed cake respectively [5].
Like many other major temperate field crops, canola is particularly susceptible to different environmental stresses, particularly heat and drought [6]. Reduced biomass and chlorophyll contents due to deterioration of chloroplast structure along with reduced production of seed oil and protein contents are usual symptoms of water stress in B. napus [7]. The plant growth and development mainly restricted by the drought stress via declining turgor pressure of plant cells which caused hindrance in the biochemical and physiological intriguing mechanisms [8]. The water deficit condition predominates with the reduced concentration of intracellular CO 2 , chlorophyll destruction, photochemical system disorder and stomatal closure [9]. The plant showed biochemical response in the form of reactive oxygen species (ROS) production under drought stress. The excess amount of ROS can damage the cell membrane by elevating the lipid peroxidation [10,11]. The plants mitigate the ROS harmful effect by the activities of enzymatic and non-enzymatic antioxidants [12]. Superoxide dismutase (SOD), peroxidase (POX), catalase (CAT) and polyphenol peroxidase (PPO) are the antioxidant enzymes which diminished the ROS concentration during drought stress [10,13]. The accumulation of osmolytes, proline, soluble sugar, soluble proteins are the non-enzymatic reaction that caused osmatic regulation under water deficient condition stress [14]. The antioxidant activities elevated under stress condition in ajowan and canola plant against ROS [15]. The plants respond to drought by enhancing osmolytes (proline), antioxidant activities (POX, PPO, SOD) in plants [16,17].
N play a pivotal role in plant tissue growth and development, being an integral part of protein, chlorophyll, nucleotides, protein, amino acid which directly influence the quality and quantity of crop production [18]. The adequate N supply is important attribute to boost up the canola productivity [19]. Thus, any fluctuation in the soil profile, texture, and moisture content at various critical stage of growth and development may decrease the N use efficiency in canola. Canola crop is very responsive to fertilizer application especially N,P and K which significantly effects the growth and yield per ha [20]. It also stimulate the leaf area (LA) development after flowering in canola [21,22].
The management of soil fertility is limiting factor for the sustainable agriculture production [20]. This challenges of soil fertility can be overcome by the optimal application of fertilizers and plant growth regulators [23]. The chemical and bio-fertilizer are very effective in improving the micro-and macro nutrient via organic compound degradation and N fixation [24]. It improves nutrient uptake and reduces the damaging effects of drought in crops. This would increase the activities of PPO, CAT and POX which ultimately improve the grain yield under drought stress in response to fertilizer and PGRs [10]. The combined application of chemical fertilizers, bio fertilizers and PGRs enhanced the accumulation of proline, sugar contents and chlorophyll contents [25]. The fertilizers along with PGRs elevated the stomatal conductance, water contents and total chlorophyll content under water limiting conditions [26]. These combined application of nutrients improve the soil fertility under water stress condition [27]. Although genetic manipulation has a promising effects to mitigate this problem and also stated mid-and long-term results, but the current demands of food and feed required some immediate response methods to address food security and hidden hunger. In spite of numerous scientific efforts, information related to biochemical, physiological and yield dynamics with respect to combined application of chemical fertilizers and PGRs is very limited for canola production. Thus, this manuscript comprised to gauge the morphological, physiological and biochemical rejoinders of canola to assimilate fertilizers management and PGRs with limited supply of water.

Experimental design and treatments
The two-year (2018-19 & 2019-2020) field experiment were conducted at experimental field of Nuclear Institute of Agriculture (NIA), TandoJam, Sindh to investigate the variability in physio-biochemical parameters and yield of canola (Surhan-2012) in response to fertilization under excessive rate of water evaporation. The source of N, P and K in the current experiment was ammonium nitrate (NH 4 NO 2 ), single super phosphate (SSP) and potassium sulphate (K 2 SO 4 ), respectively. NPK fertilizers were applied in the rhizosphere of canola growing field at the time of sowing and PGRs were applied at its flowering stage. IAA and GA 3 was thoroughly dissolved in dimethyl sulfoxide (DMSO) along with tween-20 (T-20) as surfactant. The experiment was performed in split plots as random complete block design (RCBD) in triplicate with four irrigation levels. Weather data was carefully monitored from October to April for both growing years (2018-19/2019-20). The maximum and minimum temperature decreases from October to March and then increases marginally in April for both growing seasons ( Table 1). The soil of experimental area was analyzed and found 2% organic matter (OM) along with sodium (0.08), potassium (0.11), magnesium (3.2), sulpher (16.90), zinc (1.40), boron (0.34), phosphorus (10.83) and calcium (2.5) as major nutrients. The soil mineral contents were calculated as μg/g soil ( Table 2). The detail of treatment is given in the table below.
Weeding operations was done frequently during growth and development of surhan-2012 plants in both years.

Determination of agronomic and yield related parameters
At vegetative stage, height of ten randomly selected canola plants of each treatment were measured with measuring tape (cm) from ground to the tip of flag leaf and counted number of branches per plant while number of siliqua plant -1 , siliqua length (cm) and number of seeds siliqua -1 were counted after harvesting. Similarly, recorded the data related to days to flower when 50% flowers had been appeared on plants and days to maturity at discoloration stage. Grains of canola plants were dried in the sun and recorded 1000 grain weight by using digital balance (Model-Explorer OHAUS). Biological yield was determined after harvesting from 4 rows and calculated as kg ha -1 .

Determination of quality traits
Oil content were extracted with petroleum ether using soxhlet apparatus. All the dried seed samples were coarsely ground and packed carefully into the thimble for oil extraction. The extraction was performed continuously for three cycles (90-120 min.) and oil productivity was drawn through standard formula. Two significant fatty acids, erucic acid (%) and glusinolate (μmol g -1 ) were detected through High Performance Liquid Chromatography (HLPC). Oil contents were analysed using Gradient HPLC (Shimadzu, Japan) having LC-10AT, SCTL 10A system controller, SPD-10AR UV-VIS detector at 280 nm with C18 stationary column (Shim-Pack CLC-ODS). Elution was done for 60 min with a flow rate of 1ml/min in a gradient system of two mobile phases A (H 2 O 2 : AA-94:6, pH 2.27), B (ACN100%) [28]. Moisture contents (%) was determined by the weight of water in a seed.

Determination of NPK uptake by the canola grain
For determination of P and K contents, seed samples of each treatment were dried at 70˚C for 48h. Dried and powdered grain sample (0.5g) was digested with 20mL concentrated nitric acid (HNO 3 ) by adopting method of Rathje and Jackson [29]. The samples were placed for 3 hours at room temperature. After 3 hours, samples were laid on the digestion block at 250˚C until the solution became tinted yellow in appearance. The digested solution was diluted with 50mL of distilled water and filtered with whatman No. 42 filter paper. The P contents from the digested plant samples was determined by recording optical density at 430 nm with spectrophotometer (Model-Spectronic-21) by Primson et al. [30]. The K content in gains was resolute  [31]. Nitrogen content was examined by Kjeldahl apparatus [32]. Following formulas were applied to determined NPK uptake in grains of canola as kg ha -1 . N uptake (kg ha -1 ) in grain: N (%) x grain yield (kg ha -1 ) /100 P uptake (kg ha -1 ) in grain: P (%) x grain yield (kg ha -1 ) /100 K uptake (kg ha -1 ) in grain: K (%) x grain yield (kg ha -1 ) /100 2.5. Determination of physiological and biochemical parameters 2.5.1. Chlorophyll contents. Fresh leaf samples were collected from each treatment and subjected to grinding with 80% acetone. Semi-liquid extract was filtered and centrifuged at 10000rpm for 5minutes [33]. The supernatant was then subjected to spectrophotometer (Model Analytikjena Spekol 1500 Germany) 2.5.2. Leaf water content. Leaf water content was measured by harvesting three leaves per plant from every plot after 45 days of sowing (DAS). Fresh leaf sample was weighed in gram (g) as fresh weight (FW) and let them dry at high temperature (80˚C) and reweighed as dry weight (DW). Leaf water content (LWC) was calculated by following formula.
The leaf temperature (LT) was measured at flowering stage with the help of infrared thermometer (TES-1327). The leave temperature ( 0 C) was measured by randomly selecting 3 plants of every treatment and replicate. Later, the mean LT was carefully recorded.

Stomatal conductance (g s ).
Portable photosystem (Porometer AP4, Delta-T Devices Ltd., Cambridge, U.K.) was used to measure the stomatal conductance. This data was carefully recorded 60 days after sowing (DAS). This measurement was carried out from 10:00 to 14.00 h.

Membrane stability index (MSI).
The previously reported method of Ghassemi-Golezani et al. [34] with slight modification was used to calculate membrane stability index. 0.1 g leaf samples was mixed with double distilled water (10 ml) in falcon tube and incubated at 40˚C for 30 min and electrical conductivity was measured (EC 1 ). Thereafter conductivity of these sample were assessed after placing water bath at 100˚C for 10 min (EC 2 ). The MSI was measured by the following formulas:

Determination of osmolytes
The total soluble sugar content was estimated from the dried leaves of all the replicates of respective canola treatments [35]. The standard calibration curve of pure glucose was used to determine total soluble sugars of leaves and expressed as mg/g DW. To determine proline contents in canola, leaf sample was thoroughly grinded in 3% sulfosalicylic acid. The extracted sample was filitered, mixed with glacial acetic acid and ninhydrin in a test tube with a ratio of 1:1:1. This mixture was heated at 100˚C for 60 min in a Bain Marie oven. Then reaction mixture was cooled at room temperature and the toluene used for the extraction of mixture, vortexed for 30sec. The absorbance of the upper organic phase was recorded at 520 nm. Calibration curve of pure proline was used to compare the proline content of canola leaves and expressed as mg/g FW [36].

Determination of antioxidants
Young leaves were collected from each treatment at 60 DAS and assayed the activity of polyphenol oxidase (PPO) by Kumar and Khan (1982) method [37]. The reaction mixture contains 0.1 M phosphate buffer (pH 7.8), 1 ml catechol and 5 ml enzyme extract. The reaction mixture was incubated at 25˚C for 5 min, later, the reaction was terminated by dissolving 1 ml of 3 ml NH 2 SO 4. The PPO activity was determined in the form of absorbance of resultant purpurogallin at 495 nm and expressed as Umg -1 (U = change in 0.1 absorbance min -1 , mg -1 protein. The CAT activity was determined with an interval of 20 seconds for 2 minutes at 240 nm (Ug -1 FW) according to the devised method of Singh and Sharma [38]. The POX activity was observed with an interval of 30 sec for 2 minutes at 470 nm due to guaiacol oxidation. The activity was determined from reaction solution consisted 1 ml of 1% guaiacol, 0.3 ml of enzyme extract, 2.5 ml of 50 mM potassium buffer (pH = 7.0) and 1 ml of 1% H 2 O 2 for 2 min in reaction mix [39]. The SOD activity was assessed by the estimation of volume of enzyme affected as 50% inhibition of nitroblue tetrazolium [39].

Determination of lipid peroxidation
Malondialdehyde content (mmol g -1 FW) from canola leaves was determined 60 DAS to estimate rate of lipid peroxidation [40].˜0.5 g of fresh leaves was homogenized in 5% trichloroacetic acid (5 ml), heated at 25˚C for 10 minutes and centrifuged at 1800g. The 2-thiobarbituric acid (TBA) was added in supernatant, placed at 98˚C for 10 min and cooled at room temperature. Finally, recorded the absorbance at 532 nm with spectrophotometer.

Statistical analysis
All the experimental data was recorded and subjected to analysis of variance (ANOVA) with linear models of statistics to observe statistical significant/non-significant differences among different traits of Brassica napus through computer program, Student Edition of Statistix (SWX), Version 8.1 (Analytical Software, 2005). Moreover, least significant difference (LSD) test was applied to verify the level of significance (5%) among different combination means [41].

Results
The results of canola presented in this manuscript was recorded for two consecutive years i.e. 2018-2019 and 2019-2020. The mean of all attributes have been tabulated and described in the result section (Tables 3-7).

Agronomic and yield performance of canola with fertilizers and PGRs under drought
The mean data of two consecutive years of agronomic as well as yield attributes (plant height, days to flower, number of branches per plant, number of seed per plant, biological yield per plot, 1000 seed weight and seed yield) presented a significant (p<0.01) interaction of irrigation to that of NPK and PGR ( Table 3). The plant growth was affected by the severe water stress (I 4 ), when no elicitor was provided to the canola seedlings. Plant height, leaf area and number of branches per plant decreased upto 2, 5 and 15% at maximum level of evaporation. The combination of NPK and PGRs (T 3 ) enhanced the agronomic performance under severe water deficit (I 3 = 150mm evaporation) condition by improving number of seeds/plant (1.76), biological yield/plant (5.32kg), 1000 seed weight (4.32g) and seed yield/hectare (2318kg/ha) ( Table 3). It was observed that days to flower decreased upto 1% and 8.82% with NPK and PGRs respectively under severe rate of evaporation (I 3 ). However, number of seeds per plant, biological yield and seed yield increased upto 11%, 0.4% and 28% respectively with T 3 treatment (NPK and PGRs) at maximum rate of evaporation. Grain yield was recovered with NPK (23%), PGRs (10%) and their combined treatment (17%) under least rate of evaporation (I 1 ) compared to reduction caused in non-treated canola plants (T 0 /6%). The low water supply during critical growth stage reduced the yield of canola (Table 3). Moreover, PGRs showed non-significant difference among all rates of evaporation for plant height but significantly vary for other agronomic and yield traits. The fertilizer applications significant enhanced the biological seed yield/plant (5.43kg) of canola and the highest seed yield (2450kg/ha) was recorded with NPK (T 1 ) under normal rate of evaporation (I 0 ). The T 1 and T 2 treatment presented a non-significant difference under normal irrigation condition. Moreover, data displayed a strong and significant interaction between different rate of evaporation (I) and applied elicitors (T) for all studied morphological and yield related features of canola except leaf area (Table 3).

Physiological performance of canola with fertilizers and PGRs under drought
The mean data of two-year field experiment of canola revealed a highly significant (p<0.01) response of NPK and PGRs application to chlorophyll contents under water deficit condition ( Table 4). The chlorophyll contents decline (10-67%) significantly in canola with increasing rate of evaporation (I 0 -I 3 ). However, the exogenous application of NPK (T 1 ) significantly enhanced Chl a, Chl b and total chlorophyll contents under control but combination of NPK and PGRs (T 3 ) progressively recovered the chlorophyll contents from least to severe water stress (I 1 -I 3 ). The Chl a contents decreased significantly with progression of water stress, even fortification of NPK and PGRs failed to completely mitigate adverse effects of severe rate of evaporation. A highly significant relation was observed between different rate of evaporation (I) and treatments with elicitors (T) for chlorophyll contents and also showed a recovery mechanism by promoting photosynthetic activity ( Table 4). The average of two season's data for water contents and temperature of canola leaves showed a significant (p<0.01) interaction of irrigation regimes and elicitors (Table 4). Both these traits worked antagonistically as decrease in leaf water content (LWC) ensured increased leaf temperature (LT) under severe drought (I 3 ) effect. Leaf temperature enhanced (40-6-%) in parallel to decrease in LWC (7-21%) with moderate (I 2 ) and severe rate of evaporation (I 3 ) but NPK and PGRs minimized the impact of evaporation and enhanced these features compared to their stress condition ( Table 4). The irrigation regimes I 0 and I 1 showed non-significant difference on water contents and leaf temperature and same effects were observed with NPK (T 1 ) under mild (I 1 ) and moderate (I 2 ) rate of evaporation ( Table 4).
The NPK application significantly (p<0.01) influenced the stomatal conductance of canola plant with different rate of evaporation ( Table 4). The water stress reduced the stomatal conductance (50%) under severe rate of evaporation (I 4 ) but exhibited non-significant difference under the mild (I 1 ) and moderate (I 2 ) rate of evaporation. The combined effect of PGRs and NPK (T 3 ) enhanced (1-2%) the stomatal conductance under normal (I 0 ) and mild (I 1 ) water

Qualitative performance of canola with fertilizers and PGRs under drought
The mean of two years statistical data depicted that the interactive effect of irrigation regimes and fertilizers significantly influenced the quality related traits (oil contents, oil yield, protein, glucosinolates, and moisture and erucic acid contents) of canola (Table 5). These features declined with the progression of water deficit condition. The experimental results defined the positive influence of NPK and PGRs (T 3 ) on oil yield (67-83%), oil (7-11%), protein (16-20%) and moisture (9-14%) contents while reduced glucosinolates (38-55%) and erucic acid contents (20-22%) with ongoing increasing rate of evaporation (I 1 -I 3 ) (Table 5). It is summarized that T 3 treatment is a good rehabilitation strategy to improve quality of canola followed by T 1 and T 2 to address water scarcity issues of canola.

Nitrogen, phosphorus and potassium (NPK) contents and their uptake in canola
The mean data of two consecutive years (2018-19/2019-20) showed that combined effect of fertilizer and PGRs (T 3 ) significantly (p<0.01) influence the NPK contents and their uptake in canola plant. Provision of NPK (T 1 ) individually or in combination with PGRs (T 3 ) improved the NPK contents as well as their uptake under mild (I 1 ) irrigation stress ( Table 6). The N uptake and percentage (%) accumulation has been decreased upto 22 and 5% while P uptake and % accumulation decreased upto 11 and 3% respectively, and K uptake and % accumulation halso reduced upto 17 and 14% respectively, at severe rate of evaporation (I 3 ). But plants treated with NPK and PGRs (T 3 ) enhanced overall nutrient pool (N, P, K) at all rate of evaporation (I 1 -I 3 )

Biochemical performance of canola with fertilizers and PGRs under drought
The results obtained from two-year trails illustrated that osmolytes considerably (p<0.01) influenced by the applied elicitors under different irrigation regimes (I 0 -I 3 ). The concentration of proline and total soluble sugar elevated upto 6-12% and 4-30 in response to mild to severe rate of evaporation in canola plant as an innate response mechanism. Further, the application of NPK and PGRs (T 3 ) showed a negative impact on the proline accumulation but positively enhanced soluble sugar contents with progression of evaporation rate ( Table 7). The application of NPK and PGRs (T 3 ) showed 115% increase for TSS contents at severe rate of evaporation (I 3 ). A similar trend for TSS was observed with mild to moderate rate of evaporation by the application of NPK and PGRs (T 3 ). The antioxidant enzymes and MDA activities significantly (p<0.01) influenced by the irrigation (I) and elicitors (T) in canola ( Table 7). The canola plant enhanced the activities of various antioxidants and enzymes including SOD (138-643%), CAT (182-1147%), PPO (67-304%), POX (134-1752%) and MDA (47-208%) contents in response to limited water supply (I 1 -I 3 ). The canola behavior with fertilizers did not show significant effect during normal (I 0 ) Table 6

Irrigation
Treatments N (%) P (%) K (%) N uptake (kg ha 1 ) P uptake (kg ha 1 ) K uptake (kg ha 1 ) Note: I 0 = 60mm evaporation, I 1 = 100mm evaporation, I 2 = 120mm evaporation, I 3 = 150mm evaporation, T 0 = No treatment, T 1 = NPK, T 2 = PGRs, T 3 = (T 1 +T 2 ). The alphabetical superscript in a column present significant difference among the treatments to different rate of evaporation from highest to lowest value (a = highest value). � = least significant �� significant ��� highly significant. https://doi.org/10.1371/journal.pone.0260960.t006 and mild irrigation (I 1 ) stress except lipid peroxidation and POX activity. The T 1 treatment exhibited non-significant difference on enzyme activities but the MDA content influenced significantly under all irrigation levels (I 0 -I 3 ). The treatment T 2 (PGRs) and T 3 (NPK and PGRs) elevated the activities SOD, PPO, POX and SOD but reduced the MDA contents as compared to T 0 under severe water deficit (I 4 ) condition. (Table 7). The interaction of irrigation intervals (I) and elicitors (T) significantly affect the membrane stability index (MSI) of canola (p<0.01). The water stress significantly reduced the membrane stability index (1-26%). NPK did not improve MSI while the combined application (T 3 ) significantly improve the MSI under moderate (I 2 ) and severe water stress (I 3 ) (Table 7).

Discussion
The findings of the current work highlighted the comparison of combined and individual application of PGRs and NPK to enhance canola growth, nutritional quality and yield under growing concerns of water scarcity. Particularly moderate to severe drought stress (I 2 -I 3 ) imparts drastic effect on the canola growth by inducing injuries at all growth stages. It influenced the various morphological (reduced leaf growth, leaf area, plant height, number of nodes per plant), physiological traits (chlorophyll content, leaf water contents, leaf temperature and stomatal conductance) at the onset of water scarcity (Tables 3 and 4). Reduction of plant height was recorded in canola under different irrigation regimes compared to exogenously applied PGRs and NPK (Table 3). Growth retardation due to excessive evaporation, can be related to disruption of photosynthetic machinery and decline in carbon reserves for relocation to growing parts of plant [42,43]. It seems that the decrease in plant height also interferes with leaf area. This reduction is particularly noticeable during post vegetative stage, flowering stage or abscission [43]. The grain yield reduced due to water stress as reported in pervious study [44]. Seed quality and 100 seed weight are commercial as well as economic traits, significantly compromised with ongoing scenario of water scarcity, and so is in the present experiment (Table 3). Seed filling is particularly influenced by drought stress by modulating of various metabolic activities occurring in the leaves, such as synthesis and translocation of photoassimilates, essential substrates for biosynthesis of seed storage reserves, mineral nutrients and many more functional constituents [45].
In the current experiment, drought led to reduction in chlorophyll content, and this loss could be due to some devastating effects on photosynthetic apparatus (Table 4). Stomatal conductance (g s ) was severely hampered when plants were exposed to severe rate of evaporation (I 3 ). The resistance in stomatal conductance (g s ) may be correlated to enhanced production of ABA under drought stress, which leads to stomata closure. ABA signaling mechanism tries to prevent the loss of tissue turgor by closing the stomata [46,47]. The optimal use of chemical fertilizers and PGR appreciably enhanced stomatal conductance under drought stress (Table 4). Yan et al. [48] also reported diffusional restrictions of CO 2 by stomata (52%), which directly caused a reduction of chlorophyll contents (31%) induced by drought. The fertilizer application especially urea increase the N supply at flowering and pod filling stage, delay leaf aging, enhanced chlorophyll contents and photo assimilates [49]. The total amount of chlorophyll contents increased by the availability of nitro compounds in the rhizosphere and consequently to the plants, ultimately produced more assimilates via photosynthesis which directly related to improve growth and yield [10,50]. Growth regulators and chemical fertilizers were significantly effective in mitigating the drastic effects of drought by maintaining the water efficiency of canola plants and augmenting the accretion of osmolytes. Accumulation of osmolytes may also favors the improvement of photosynthetic and gas exchange attributes [16]. The observed increase in yield of canola using NPK and PGRs under water limitation may be attributed to enhanced activities of CAT, SOD, POX and PPO [51]. Moreover, combined application of chemical fertilizer, PGRs and vermicompost particularly enhanced the accretion of secondary metabolites such as proline and sugar content and also chlorophyll synthesis [52].
Canola oil is the commercial commodity, while its content, profile and composition are affected by drought stress as reported in the current work (Table 5). Seed oil stems mostly from photosynthesis and green silique walls, later carbon is routed through different metabolic pathways into triacylglycerol occurring in the chloroplast, cytosol, and endoplasmic reticula [53]. The current experiment suggests 3-11% decrease in oil content and oil yield in the B. napus when the plants were exposed to irrigation stress, but Aslam et al. [54] reported a mere reduction of 3.2%.
The NPK treatment increases the flow of nutrient to the aerial part of plants and reduced the impact of water stress (Table 6). Due to water scarcity and loss of ionic balance, nutrients remained bond to the soil particles that are critical for the normal growth and development [6,55]. Same findings have been reported in the present work in the form of decreased biological yield (40%) with non-availability of nutrients (5-10%). Particularly, N, P and K contents increased with foliarly applied PGRs followed by chemical fertilizers in canola. This might be accredited to the role of K in biochemical pathways in plants. Potassium has a positive effect on metabolic processes of nucleic acids and proteins [56]. Phosphorus as a constituent of cell nuclei is essential for cell division and development of meristematic tissue of cotton. Further, P has a well-known impact in photosynthesis as well as synthesis of nucleic acids, proteins, lipids and other essential compounds [57]. The percentage of NPK uptake enhanced with the combined application of NPK and PGRs ( Table 5). The bio fertilizers improve the soil textures and bacterial colonization with the modification of physio-chemical properties of rhizosphere. On the other hand, the PGRs increased the photo assimilates translocation to sink (root tissue) and also improved the nutrient uptake and absorption power under water deficit and adverse environment condition. The PGRs can enhanced the activity of some vital N, P and K metabolizing enzymes in plant which enhanced NPK contents under different irrigation regimes [50,58].
All plants have been equipped with innate antioxidant enzyme mechanisms for the detoxification of reactive oxygen species. CAT decomposes H 2 O 2 into water and molecular O 2 . POX converts H 2 O 2 by oxidizing co-substrates such as phenolic compounds and/or antioxidants and PPO in turn oxidizes phenols to chinone [59].
Membrane lipid peroxidation is a frequently used indictor to test the degree of plant sensitivity to oxidative damage caused by stress [60]. Our data revealed sensitivity of canola towards drought by elevation MDA content compared to F 2 , F 3 treatments (Table 7). A lower level of lipid peroxidation presented high membrane stability (Table 7). It seems that the cell membrane integrity was maintained with chemical fertilizers and PGRs against the oxidative stress induced by water stress (Table 7). Mamnabi et al. [10] also suggested an amplification of antioxidant enzymes, total soluble sugars, photosynthetic attributes, leaf water content, membrane stability index and stomatal conductance but decreased the leaf temperature under different irrigation regimes.
The affirmative role of fertilizers and PGRs increase the antioxidant enzyme activities under water deficient condition in canola [61]. The PGRs treated plant showed more activities of antioxidant enzymes like POD, CAT and PPO as compared to untreated plants under moderate and severe water stress. The highest antioxidant activities were observed in T 3 treatments under severe stress as compared to T 1 and T 2 ( Table 7). This supremacy was attained by the additive effect of fertilizer and PGRS on canola plants.

Conclusion
Canola is an emerging and unique oil seed crop among the other oil producing plants due to low erucic acid and glucosinolate contents. Here, we initially illustrate how NPK and PGRs, either individually or in combination, impact canola growth, photosynthetic and antioxidant activities and later seed yield and quality, and also attempt to explain its interaction to water scarcity for addressing these vital challenges. From the outcomes of current study, it appears rational to recommend chemical fertilizers (NPK) and PGRs (IAA and GA 3 ), that brought about better impact on canola seed yield, seed protein content, oil, oil composition with low glucosinolate and erucic acid contents, even under severe rate of evaporation (150mm). The harmful effects of stress were minimized considerably by the combined application of fertilizers and PGRs, thus improved growth attributes, chlorophyll content, MSI, stomatal conductance, antioxidants, osmoprotectants, grain yield and importantly, leading to a reduction in lipid peroxidation, particularly under moderate and severe rate of evaporation. These supremacies were attained by additive effects of NPK and PGPR, reducing the impact of drought. Such models can improve the probability of forecasting canola aptitude in challenging climates with an immediate response, but will also broadly help to select traits that can be further exploited through gene mining to produce sustainable and climate-resilient canola genotypes with considerable yield under high rate of water evaporation.